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THINGS HAYING ALL EDGIS SUBSONIC 
By Harold J. Walker and Mary B. Ballantyne 


- SUMMARY 

A method is presented for calculating the pressure distribution and. 
damping in steady pitch at supersonic Mach numbers of thin, flat, swept— 
back wings having all edges straight and subsonic. Although it is adapt- 
able to wings with negative rake at the tips, the method is applied only 
to wings with streamwise tips. 

The method consists of two steps: first, the calculation of a basic 
pressure distribution, which is identical to that existing on an infinite 
triangular plan form having leading edges that coincide with those of the 
swept-back wing; and second, the correction of this basic distribution to 
account for the effects of the subsonic trailing edges and tips. In cal- 
culating the various corrections, use is made of the principle of the 
superposition of conical flows. The derivative for the damping in pitch 
is calculated in a similar manner. 

In applying the method of analysis to a typical configuration, it 
was found that several of the corrections .were small, and, from a practi- 
cal standpoint, could be dropped. By dropping these terms the method is 
shortened to the extent that it closely parallels previously published 
methods for calculating the lift, pitching moment, and damping in roll of 
swept-back wings with subsonic edges. 

A substantial reduction of the pressure in the vicinity of a subsonic 
tip was disclosed in the analysis. This effect was also found earlier for 
the cases of steady lift and steady roll. 

The method is based upon the usual assumptions and limitations of 
the linearized potential theory for supersonic flow. • 


INTRODUCTION 

Methods have been developed previously for the calculation of the 
pressure distribution and the damping in steady pitch at supersonic Mach 
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numbers far thin, flat wings of various plan forms, including triangular 
wings (reference l), rectangular wings (reference 2), wings with super- 
sonic leading edges hut somewhat arbitrary plan form (reference 3)> and 
swept-back wings having all supersonic edges or a combination of subsonic 
leading edges and supersonic trailing edges and tips (references h and 5). 1 
An approximate analysis of the damping in pitch of a l imi ted class of 
swept-back wings having all edges subsonic and straight is reported in 
reference 6. In the present report, a method, which can be applied more 
generally than that of reference 6, is presented for the analysis of swept- 
back wings with straight subsonic edges. 

The present method of analysis consists essentially in the super- 
position of various coni cal and qrasi— conical 2 fields of pressure in such 
a mann er that the particular boundary conditions for a swept-back wing in 
steady pitch are fulfilled. In detail, it closely follows the analyses 
given in references 7 &nd 8 far the case of steady lift and in reference 9 
far the case of steady roll. The analyses far the three cases, in fact, 
differ only in application to the different sets of boundary conditions 
corresponding to lift, roll, and pitch. 

The present analysis and those of references 7 and 9 are limited to 
wings having zero or negative rake at the tips; however, wings having 
positive rake can be treated through an adaptation of the method given 
in reference 8. Although the analysis for the type of configuration in 
which the Mach lines from the root of the trailing edge intersect the 
leading edge ip not complete, the results given for this case are suffi- 
ciently accurate for many practical purposes. The scope of the method 
in general is limited to the usual idealizations and assumptions of lin- 
earized potential theory. 

To illustrate the application of the method, calculations of the 
damping derivative and of the pressure distribution along several chord- 
wise and spanwise sections of the swept-back wing Bhown in figure 1(a) 
are included. 


DOTATION 


a slope of any ray through origin divided by slope of Mach 

lines (fig. 1(b)) (pj) 

1 The terms "subsonic edge" and "supersonic edge" refer to edges having 
normal components of flow which are subsonic and supersonic, 
respectively. 

The term "quasi— conical pressure" designates a pressure which is 
distributed linearly along rays passing through a fixed point. 
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Co 


the ray a intersecting the trailing edge at the point 
from which the Mach line passing through the tip of the 
leading edge emanates (fig. 1(b)) 


a Z ~ 




1 - 


Pb 

Conmifc 


-(1-m) 


•the ray a intersecting the trailing edge at the point 
from which the Mach line passing through point P(x,y) 

r py+c D — x 

emanates (fig. 1(h)) ! on trailing edge, a 0 =mj- ; 

Py+Com-t-x 


on tip. 


a o= 


1EL 


x+P ( y-s ) 




slope of ray through the trailing— edge tip divided hy 
slope of Mach lines (fig. 1(h)) ^ ^g£L ) 

wing Bpan 


“t 


pitching-moment coefficient 



derivative for the damping in pitch 

(In the present report, with no subscript represents 

the damping moment about the c enter— of —gravity axis 
resulting from the loading due to steady pitch about the 
y axis instead of the center— of— gravity axis. See 
equation (2).) 



root chord of wing 


* 5 

m 

“a 


auxiliary functions used in analysis 

slope of leading edge divided by slope of Mach lines 
(fig. 1(b)) (p cot A) 

slope of ray passing through point x A ,y A and tip of lead- 
ing edge divided by slope of Mach lines (fig. 1(b)) 
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slope of ray passing through root of trailing edge and 
tip of leading edge divided "by slope of Mach lines 

■ <«•■««> (pg;) 

slope of tip divided "by slope of Mach lines 

slope of trailing edge divided by slope of Mach lines 
(fig. 1(b)) 

pitching moment about the y axis 

free-stream Mach number 

pressure difference between the upper and lower wing sur- 
faces 

Ap 

pF 

correction to the basic pressure coefficient 

basic pressure coefficient at the root of the trailing 
edge for a wing at constant angle of attack 
* 

basic pressure coefficient at the root of the trailing 
edge for a wing in steady pitching motion 

steady pitching velocity, radians per unit of time 

(positive, if direction of rotation is that of increas- 
ing angle of attack) 


pressure coefficient 


( 


1 — 2m 2 xr » m 2 -ci 

: 2 K + : a E 

1 — m 1 — m 

semispan of wing 

area of wing 

slope of any ray passing through the point, x A ,y A 
divided by slope of Mach lines' 


s ( P — ) 
\ y-y A / 


slope of ray through trailing— edge apex divided by slope 


of Mach lines 


(o^) 
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x,y,z 
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X 
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• 6 * 


Z 

qc 0 

2V 


P 

8 

P 

X 

A 



free-stream Telocity- 

Cartesian coordinates in the stream direction, across the 
stream, and in the vertical direction, respectively (fig. 
1(D)) 

streamwise distance from y axis to the center of pressure 
of an element of wing area 

streamwise coordinate of point where any ray t Q intersects 
tip or leading edge (for the tip, x s = + c Q ; far the 

c~m . Do 

leading edge, x s = + c ) 

t Q — m o 

streamwise distance from y axis to center of gravity of 
airplane configuration 

force in the vertical direction 

dimensionless parameter representing steady pitching velocity 

'/m q 2 -i 

constant factor used in equations (5) and (6) 
free-stream mass density 

argument of inverse cosine terms (see text) 
angle of sweep of leading edge 

Elliptic Integrals 

complete, elliptic integral of first kind with modulus J 1-m 2 
complete elliptic integral of second kind with modulus 1-m 2 
complete elliptic integral of first kind with modulus k Q 
complete elliptic integral of second kinfl with modulus k D 

incomplete elliptic integral of first kind with modulus k Q , 
amplitude cp Q 


F(<P 0 ,k 0 ) 
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E(<p 0 ,k 0 ) incomplete elliptic integral of second kind vith modulus k Q> 
amplitude <p Q 


sin" 1 

71 -mt 2 

complete elliptic integral of first kind with modulus k 
complete elliptic integral of second kind with modulus k 

/ m-ap) (1-m) 

2m(l+a 0 ) 


incomplete elliptic integral of first kind with modulus k* , 
amplitude cp 

incomplete elliptic integral of second kind with modulus k*,- 
amplitude 9 


1 “ J 

k« V 2 


incomplete elliptic integral of first kind with modulus k* , 
amplitude ijr 

incomplete elliptic integral of second kind with modulus k*, 
amplitude ijr 


sin -1 i 

k * V Py+x 2m 

[E(k) -K(k)] P(q>,k«) + K(k) E(q>,k’) 
[E(k) — K(k)] F(t,k») + K(k) E(i»f,k*) 


Subscripts 


2 


terms related to the conical pressures (except as noted) 
terms related to the quasi-conical pressures (except as noted) 
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A 


L 

P 

R 


a 


terms evaluated at the point A(x A ,y^) 

on the trailing edge , xa = j 

1 m^-a A 

terms corresponding to steady lift 

terms corresponding to steady pitching velocity 

termB corresponding to steady rolling velocity 

terms referring to center— of-gravity axis of airplane 
configuration 

designates "per unit angle of attack* 1 



on the wing boundary 


_ m t c o a oq the tip 
3(mfc-a) 


Superscripts 

terms related to symmetrical canceling sectors in wake 
terms related to oblique canceling sectors in wake 
terms related to canceling sectors outboard of tip 

METHOD 0? ANALYSIS 


The method for calculating the pressure? and the damping derivative 
for flat swept-back wings in steady pitch is developed in a manner simi- 
lar to that in reference 7 for wings in steady lift anfl in reference 9 
for wings in steady roll. The present analysis, therefore, will be 
shortened somewhat by referring frequently to these two reports. Tha 
analysis is restricted to wings having straight subsonic edges and, 
although it can be adapted easily to wings with negative rake 4 at the 
tips, is applied only to cases of zero rake (i.e., streamwise tips). 

A thorough treatment of those cases in which the Mach lines originating 
at the trailing edge intersect the leading edge is not included; however, 
results which are sufficiently accurate for most practical cases of this 
type can be obtained from the method as outlined. 

3 

Throughout the report the terms "pressure" and "pressure coefficient" 
will be used synonymously. 

The rake at the tip is negative if the tip slopes inwardly from the 
leading edge. 
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In the coordinate system chosen for the analysis, the origin is 
located at the apex of the leading edge and the x axis coincides with 
the root-chord line of the wing. (See fig. 1.) The y and z axes 
extend spanwise in the plane of the wing and perpendicular to the plane 
of the wing, respectively. 

The analysis is facilitated hy assuming initially that the pitching 
axis coincides with the y axis. Subsequently, the pressure distribution 
due to pitching about an axis passing through the center of gravity of the 
airplane configuration (assumed to lie an the x axis at a distance 
x downstream from the y axis) can be obtained by superposing on the 
pitching motion about the y axis a vertical translational motion of 
velocity qx c< g # o The pressure distribution due to this translational 
motion corresponds to that of a wing at a c ons tant an g le of attack of 
-qx c .g./V. (See reference 4.) Thus, 



where Pp and have been corrected for the effects of canceling 

the excess basic pressure in the wake and outboard of the tips. The cor- 
responding transf orrnati on for the damping derivative is 


(si “ (O 


c.g 
in which 


are the damping derivative and lift 


coefficient corresponding to steady pitching motion about the y axis, 
and and Cj ^ are the moment coefficient and lift coefficient per 

unit angle of attack. In order to simplify the analysis in the present 

are calculated as a single 


report, the terms and ^Cj,^ 


term C m such that 

q 


(c m '\ = C m - 2 Xc -&- C™ -2 Or 

V “q °b "k V c o / ^ 

c.g. 


( 2 ) 
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It should he observed that the term then represents the damping 

moment about the center-of-gravity axis resulting from the loading due 
to steady pitch about the y axis, and therefore does not, by itself, 
constitute the c cample te damping derivative (except in the spec ial case 
in which the center of gravity lies at the leading-edge apex) . The 
additional terms (Pj^, C-^, and C^) required for other locations 


of the center of gravity can be calculated by the method of references 
7 and 8. 


The order of the analysis consists in the derivation of the pressure 
distribution, followed by the calculation of the .damping derivative. 


Pressure Distribution 


Following the procedure given in references 7, 8, and 9, the swept— 
back wing is considered initially to be an integral part of an infinite 
triangular plan form, the leading edges of which coincide with those of 
the swept-back wing. Then it is possible to calculate a basic pressure 
distribution over the swept-back wing by means of the simple expression 
for the tria n gula r wing. The excess basic pressures introduced in the 
wake and outboard of the tips of the swept-back wing in the first step 
are subsequently canceled by superposing over those regions a series of 
conical and quasi— conical pressure fields. At the same time, the Kutta 
requirement of zero pressure along the trailing edges and tips is ful- 
filled. If the trailing edges and tips are supersonic, these exterior 
pressures may be neglected since their influence does not extend onto 
the plane of the wing. For wings with subsonic edges, several primary 
and secondary corrections representing the effects of canceling the 
excess pressure in the wake and outboard of the tips must be added to the 
basic pressures. 

Basic pressure distribution .— The basic pressure distribution far 
the trian g ula r wing in steady pitch is given in reference 4 as 


P P = 



8 x 2m 2 — a 2 
P® C ° /m 2 -a 2 


B = 


l-Sa 2 

1-m 2 


K + 


m 2 

1-m 2 


E 


(3) 


The variation of pressure is seen to be quasi conical (i.e., the pressure 
varies linearly in the x direction along a given ray a) and to be 
dependent principally upon the sweep of the leading edge. Equation (3) 
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also can "be used directly to calculate the pressure distribution of swept— 
back configurations, except in the vicinity of a subsonic bo undar y. For 
the configuration considered herein, the regions in which the above 
expression is inadequate are labeled I, H, anfl hi in figure 1(a). 

Region I is influenced by the excess pressure in the wake, region H by 
the excess pressure outboard of the tip, and region HI by the pressure 
both in the wake and outboard of the tip. 

Primary corrections due to cancellation of basic pressure in wake.— 

The corrective terms, which in regions I, II, and HI correspond to the 
induced effects caused by the cancellation of the excess basic pressure, 
are designated "primary corrections" in the following analysis. They are 
derived by superposing along the subsonic boundaries a series of sectors . 
of conical and quasi— conical pressure which, when integrated, completely 
cancel the excess pressure. The function for each sector of canceling 
pressure, in order to fulfill the boundary conditions of the pressure to 
be canceled, must-; 

1. Represent a field of pressure consisting of conical and quasi— 

conical portions which conform with the pressure field given 
• by equation (3) 

2. Have a downwash flow field which does not, in general, extend 

onto the plane of the wing, in order that the flatness of the 
wing be maintained 

3. Satisfy the linearized equation for potential flow 

A function, which fulfills the above requirements with regard to the 
cancellation of the basic pressure in the wake, iB composed of two prin- 
cipal parts: the first, represented by a single symmetrical sector of 
pressure, and the second, represented by a series of oblique sectors of 
infinitesimal pressure. The first cancels the relatively large field of 
pressure determined by the basic pressure at the apex- of the trailing 
edge (fig. 2(a)), and the second cancels the remaining smaller portion 
of the pressure (fig. 2(b)). 

The first principal part of the function is composed of a conical 
and a quasi— conical component, and can be derived from the symmetrical 
function previously utilized in reference 7 for the lifting case. Thus 
the function 5 


5 It is understood that the real part of F(<P 0 ,k 0 ) applies. Far values 
of sin <? greater than one, the real part of cp is equal to rt/2, 

O c O 

and the real part of F(cp 0 ,k 0 ) is equal to K Q . 
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£P t = - 


Ko 




where 


k Q = V 1-n^ 2 



and, where P T is the lifting pressure at the apex of the trailing 
J -'o 

edge, can he used directly in the pitching case to cancel a field of 
constant pressure, equal in magnitude to the pressure (due to pitch) 
at the trailing-edge apex, hy rewriting it in the form 


APp*’ - - T 2 *<V* 0 ) 


w 




1 - in t' 
<Po = sin 


■i /EZ 

J l-m* 2 


where Pp Q is given hy equation (3) for a=0 and x=c Q , that is, 

p =( 2£o\ 1^ 

P o V 27 / P® 

The field of pressure canceled hy means of equation (4) is conical 
(fig. 2(a)), since the canceling pressure is constant along the rays t Q 
originating at the apex of the trailing edge. 


A symmetrical function, which cancels the remaining quasi-oonical 
portion of the pressure (fig. 2(a)), can he derived hy the method of 
superposing infinitesimal conical, pressure fields, as outlined in 
reference 9* The derivation consists in the formation of a function for 
a quasi— conical pressure field hy integrating the effects at a point 
P(x,y) on the wing of a series of infinitesimal symmetrical conical sec- 
tors represented hy equation (4) . The sectors are distributed in a pyram- 
idal arrangement with the apexes located along the ray to = 0. If | is 
the streamvd.se coordinate of the apex of each sector measured from the apex 
of the trailing edge, and the limit g 0 designates the apex of the rear- 
most sector containing the point P(x,y) within its Mach cone, the 
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derivation leads to the expression 


= r.p. 83 . / 




= r.p. 8 a 


in which 


" „ , Itol -i / m t 2 ( 1 ~^o 2 ) 

_ ; F( c P 0 J k 0 ) j^l sin y to 2 (l-Tttfc 2 ) 


i x-Cq-I 

and 83 . is a constant. As desired, the real part of equation (5) is 
directly proportional to the distance (x-c Q ) and varies linearly with 
t Q between the ray t Q = 0 and the trailing edges of the wing as shown 
in figure 2(a) . The expression, however, reduces to the value 
8 1 (x~c 0 ) (Kq-jc/2) / k q3 instead of zero, along the trailing edges. Hence, 
in order to satisfy the boundary conditions along the trailing edges, it 
is necessary to subtract from equation ( 5 ) 811 additional, function which 
has the value &! (x-c Q ) (Eq-it / 2 ) in the region of the wake. This addi- 
tional function can be derived from the two expressions obtained by super- 
posing the infinitesimal conical sectors represented by equation (4) along 
the rays t Q = mt and t Q = -mtj respectively, in a manner similar to 
that for deriving equation (5) . In thiB manner, the two functions 

f 2 =r.p. 82 (x-c 0 )[mt(l-^fct 0 )F(cp 0 ,k 0 )+(t 0 -m t )E(cp 0 ,k 0 ) - | t °^j J (l--fc 0 2 )(t 0 2 -P^ 2 )] 

and 

f 3 =r -* p * ^(*-°o) [^(i+mttoJP^o^o Hto+^Efoo^)* j t °^| y(l-t 0 g )(t 0 g -mt g ) J 

are determined. Since only the real parts of these two functions are to be 
considered, the factors which determine the signs of the radical terms may . 
be replaced by t Q / jt 0 1 . The two functions are now added giving the 
resultant function 

^ = r.p. 8 4 (x-c q ) [F(<P 0 ,k 0 ) - E(q> 0 ,k 0 ) ] 


s The abbreviation r.p. indicates real part. 


(6) 
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which, along the trailing edges reduces to the value 8^ (x-o Q ) (Kq-Eq ) . 

If 8 4 is set equal to &! (E Q — n/2 ) /E q (E 0 '-E q ) , then equation (6) 

becomes the additional function required to satisfy the boundary condi- 
tions along the trailing edges. 7 Except for the determination of the con- 
stant 8i, the expression obtained by combining equations (5) and (6) 
describes the desired field of quasi-conical canceling pressure. This 
function must have the same boundary value along the x axi s as that 
given by equation (3) with x replaced by x-c 0 , that is, 

6 n = f = p x -°o 

1 Ko 2 \ 2T / pp Co “ ^ P o °o 


Thus, 8i must be equal to 2Pp o K 0 /itc 0 , and the corrective tern corre- 
sponding to the quasi-conical portion of the pressure becomes 



gp-Wg) 


[F(<Po,ko) ~ E(tpoj^o)] ^ 


(7) 


Equations (4) and (7) together comprise the first part of the over- 
all function required to cancel the excess basic pressure in the wake. 

The field of pressure represented by these equations is shown in figure 3, 
in which the regions of overlapping (or induced) pressure on the wing are 
indicated by the dotted lines. 

The second part of the complete pressure-canceling function for the 
basic pressure in the wake (consisting of a series of oblique sectors) 
has been developed in reference 9 , but in the form, appropriate to a condi- 
tion of steady roll. The expression for a single oblique sector of cancel- 
ing pressure for the right half of a rolling wing is 


7 It is interesting to note that, if the difference between the func- 
tions fi and f 2 is taken, the function 

f 5 = r.p. |s(cp 0 ,k 0 ) - m/ F(cp 0 ,k 0 ) - ^ /(l-t 0 2 )(t 0 2 ^a,. 2 ) J 

which is directly proportional to y, is obtained. This function 
and f 4 have been derived by other methods in reference' 10, and are 
utilized in reference 6 to calculate approximate values for the damp ing 
in roll a n d pitch for swept-back wings with subsonic trailing edges. 
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iia/i-x.. 2 ) aa 

t-m^ 1-a / 

= d(AP Ei ”) +t.(& Ss ") 

(l-a)(t-mt) - (mt-a)(l-t) 


( 8 ) 


where 


X*« = 


(l-m^.)(t-a) 

and where Pp^ is the pressure in steady roll at a point A(x A ,y A ) on 
the trailing edge defined "by 

_ ^fcQp 


Za = 


mt-a 


T = m t°o a 

A p (hlj.— a ) 

Equation (8) is seen to he made up of two components, the first 


conical., and the second 


Ov) 




quasi conical [due to the term (y— y^J/y^] . 


Equation (8) defines an oblique sector of infinitesimal pressure to 
he superposed on the basic pressure in the wake. As shown in figure 1(h), 
each sector has an apes at the point A(x A ,y A ) an the trailing edge, a 
base at an infinite distance downstream, and is hounded on one side by the 
ray a wnfl on the other by the trailing edge. This expression can be 
used to cancel pressures due either to roll or pitch; however, for the 
case of pitch it is necessary to rewrite the equation in terms of x and 
x A rather than y and y^. Thus, by substituting 


x ~ x a _ ^ y ?A 
x A t y A 


and by replacing P R with Pp, equation -(8) is transformed to 

) co B - X-aa-i rst) =» =£ (co=- X. 

\ da / « \ da / x A mfc-a \ 


aUay) - -I 


t _ 


t-mfc 1-a 


da 


) ( 9 ) 


. a(a? Pi ") + a(app 2 ") 
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The field of pressure represented "by equation ( 9 ) is illustrated in fig- 
ure 4, in which the lower portion of the pressure corresponds to the coni- 
cal component and the upper portion to the quasi— conical component. The 
region of the pressure 'which overlaps the wing is shown "by the dotted lines . 


The total effect upon the pressure at a point P(x,y) on the right 
wing of canceling the field of "basic pressure in the wake in excess of the 
pressure canceled by equations ( k ) and ( 7 ) is found by integrating equa- 
tion ( 9 ) between the l imi ts a=0 and a=a Q . The upper limi t a Q repre- 
sents the rearmost sector, the Mach cone of which passes through the point 
P(x,y) (fig. 1(b)), and is defined in reference 7 as 


0y+c o -z 

a o " “t p y+mtCo _ x 


( 10 ) 


Also, for located on the trailing edge, it can be shown that 

*-*A Vl = m t Z ~^ 7 _ -l 

Z A m tf a m t°o 


Hence, for the basic pressure in region I, the second part of the primary 
correction becomes 



where 

a?P A = /qc Q \ &iifc f 2m 4 -3rn g a 2 + a 3 !^ 1 
da '2V/ PE l(m t -a) 2 (m 2 -a 2 ) 3 / 2 J 


A graphical method of integration is recommended far solving equation (ll). 
The separate conical and quasi— conical components will be retained through- 
out the remainder of the a n alysis in order to show their relative magni- 
tudes when applying the method to a typical configuration. 

Equations (4), (7), and (ll) therefore comprise the complete set of 
primary corrections to be added to the basic pressure in region I adjacent 
to the subsonic trailing edge. They are only partial corrections far 
region III, which is also affected by the cancellation of pressure out- 
board of the tips and by the cancellation of certain secondary pressures 
to be discussed later. 
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Primary corrections due, to 'basic— pressure cancellation outboard of 
tins In general, the excess basic pressure outboard of any subsonic 
boundary (except the leading edge or a tip with positive rake) may be 
canceled either wholly or in part by means of equation (9) . This expres- 
sion can be readily modified to conform with .the boundary conditions of 
the excess basic pressure outboard of the tips. The section of this 
excess pressure at the tip is shown in figure 2(b). With reference to 
figure 5 and -not ing that the wing lies in the negative range of t of 
each sector of infinitesimal canceling pressure to be superposed outboard 
of the tips, it is seen that to apply equation (9) to the tip it is neces- 
sary only to substitute -mg, — t, and -a for mj., t, and, a, 
respectively. Thus, for a right-hand tip with zero rake (mg=0). 


afapp 1 ") cos lx '" to - 


i Lb - 1 

* \ da / x A a L 


X*** - 


t-a 

t(l-ta) 


a/ i_X* * * 2 J 


da 


) ( 12 ) 


= d + d^SPp a ***^ 


where X* 1 reduces to 


and where 


Xttt = a+t+2at 
t-a 


x. 

A a 


_ f qc 0 \ Suit- l~ 2m 4 — 3m 2 a 2 +a 3 mt ~ [ 
da 2V/ PR L (m^-a) 2 (m 2 -a 2 ) s / 2 -l 


The sector of infinitesimal canceling pressure represented by equation (12) 
is shown in figure 5. In this analysis only the tip with zero rake is 
considered; however, the tip with negative rake can be treated in an analo- 
gous mann er - The correction for the tip with positive rake along which 
the pressure is infinite, although not considered herein, can be calculated 
by adapting to the present case the method given in reference 8 for the 
case of steady lift. 

The total effect on the pressure at a point P(x,y) of canceling 
all the pressure outboard of the tip, as shown in figure Is found 
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"by integrating equation (12) between the limits a=a 0 and a=m, where 
a 0 for the tip is given in reference 7 as 

a = Pb 
° x+p(y-s) 

It is noted that the term dPp /da becomes infinite at the upper limit 

of integration m; therefore , A the integration should be made according 
to the method of reference 7* This method leads to the expression 


AP t 


m 

' \ 

[■ P ?A t ( 

cos - 1 X' * * J da - 

& o 



[ X-X A t / -i 

L * A a ** 


*tr-€l 

t(l+a) 


J i-jc *»* 2 ^ 


da 


which, after substituting the expression for Pp from equation (3)> 
reduces to A 


zautt. JL / a o( s ~ r) r 

P \ 2 vj -/ s L J c o 


(Py+x)(2m g -a 2 ) 


■ da + 


(Py-ax) </ (m 2 -a 2 ) (a-a D ) ( 14a) 


/ m 

X A ^o 
c^ a_(l-H 


(Py+x) (2m 2 -a 2 ) 


c o a o( 1+a ) (Py-ax)«/ (m 2 ra 2 ) (a-a Q ) (l+a) 


da 


] 


>(13) 


= AP Pi '»» + APp 2 »»* 


J 


The first integral represents the conical portion of the canceling pres- 
sure, and the second integral, the qua si-coni cal portion. Equation (13) 
has been integrated in terms of elliptic functions in Appendix A. This 
equation therefore becomes the primary correction for the basic pressure 
at points contained within the Mach lines from the tip of the leading 
edge (region II, and region III in part). 

In Appendix A, it is interesting to note that along the Mach line 
originating at the tip of the leading edge (corresponding to a Q =m) the 
conical part of the primary correction for the tip does not reduce to 
zero. This correction therefore represents an abrupt drop in pressure 
which, as is shown in the following illustrative example > is sufficiently 
large to cancel nearly all the basic pressure between the Mach line and 
the tip. Similar effects were noted in the cases of steady lift and 
steady roll. 



18 


MCA TU 2197 


Second a ry corrections .— Thus far, "by means of equations ( 4 ), (7), 

(11), and (a 5), all the excess "basic pressure in the wake and outboard of 
the tips has been eliminated. However , since the canceling sectors used 
‘ in the process are alsq infinite in extent, still further, although much 
smaller, excess pressures are again introduced in certain regions in the 
wake and outboard of the tips. For example, as shown in figure 6 for the 
right half of the wing, excess secondary pressure will be added to the 
shaded region outboard of the tip and, in some cases, forward of the 
leading edge as a result of superposing the symmetrical and oblique can- 
celing sectors along the trailing edges. In a like nanner, excess second- 
ary pressure still remains in small regions downstream from the trailing 
edge after canceling the basic pressure outboard of the tips, as shown 
for the left half of the wing in figure 6 (the secondary pressure in 
extreme cases overlapping the opposite wing panel, as shown) . A rigorous 
mathematical analysis would require the cancellation of these secondary 
pressures and, in turn, the successively smaller excess pressure introduced 
by canceling the secondary pressures, and so on. The details of canceling 
these secondary and lesser pressures by numerical methods are discussed in 
references 7 and 8; but, since the procedure is rather complicated and 
leads to only minor improvement of the final value, approximate methods 
will be utilized here. a Thus, if the secondary pressure is neglected, 
the pressure along the portion of the tip and trailing edge affected "by 
the secondary pressure will not have been reduced to zero. The magnitude 
of this error becomes evident after calculating the chordwise pressure 
distributions along sections near the tip. The secondary correction for 
the error can then be easily estimated using, as guides, the trends of 
the primary corrections. As shown in figure 1 , the regions in which the 
secondary corrections apply are located between the wing boundaries and 
the Mach lines reflected from the points where the Mach line from the tip 
of the leading edge intersects the trailing edge, and where the Mach line 
from the apex of the trailing edge intersects the tip. The excess second- 
ary pressure adjacent to the leading edge is neglected in cases in which 
the Mach lines from the trailing— edge apex intersect the leading edge. 

This approximate procedure is sufficiently accurate from a practical 
standpoint, as will be shown in the following illustrative case. The 
effect on the pressure distribution of secondary pressure existing upstream 
from the leading edge is not treated here, but is believed, on the basis 
of the results given in reference 8 far the lifting case, to be small 
enough to be neglected. 

Violations of dcwnwash boundary conditions .— Although the equations 
for the sectors of canceling pressure in general do satisfy the boundary 
requirements far the pressure, they do not in every case comply with the 
condition that the dcwnwash flow on the wing be zero in order that the 
wing be flat. As discussed in reference the terms cos -1 X 8 * * 11 and 

8 It should be noted also that a detailed analysis of the small sec ond ary 

corrections would not be fully justified in view of the possible signi- 

ficant effects of viscosity, which are not considered in the present 

analysis. ' 
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y l-X * 12 in the expression for a Bingle oblique canceling sector for the 
wake (equation ( 9 )) have a real part, which corresponds to pressure, and 
an imaginary part, which corresponds to downwash. Far these oblique 
sectors, the real parts are zero only in the range — 1 < t < a, and the 
imaginary parts are zero only in the range < t < 1. It is observed 
that some of the rays (in the negative range of t) from the sectors 
near the apex of the trailing edge of one wing panel (say the right panel 
in fig. l) will pass over the opposite wing panel, thus introducing same 
downwash flow on that panel. The same effects occur on both panels, but 
the regions of the wing and the amount of the downwash flow involved, in 
general, are small. Hence, some inaccuracies in the cancellation proce- 
dure are unavoidable, but they are believed to be insignificant in the 
final result. 

Illustrative application.— To illustrate the application of the fore- 
going analysis, the pressure distributions along chardwise sections ArA 
and B-B and along spanwise sections C-C, D-D, and E-E of the wing in 
figure 1 have been calculated and the results plotted in figure 7. In 
this example, the y axis is assumed to be the axis of pitch. For camr- 
parison of their magnitudes, the conical and quasi— conical terms in the 
trailing-edge and tip corrections are shown individually in parts (a) and 
(b), where it is observed that, in general, the quasi— conical terms cure 
small compared to the conical terms. Both for thiB reason and because they 
involve a considerable portion of the computing time, the quasi— conical 
components probably can be dropped in most practical cases. The procedure 
for calculating the pressure distribution then becomes essentially the 
same as that outlined in references 7 and. 8 far steady lift. 

A close approximation to the pressure distribution in steady pitch, 
therefore, may be calculated by the following steps: 

1. Calculation of the basic pressure distribution for the entire 

wing by means of equation ( 3 ) 

2. Correction of the basic pressure between the subsonic trailing 

edge and the Mach line from the apex of the trailing edge 
(regions I and III) by means of equations (4), ( 7 ), and (-11) 

3 . Correction of the basic pressure between the subsonic tip and 

the Mach line from the tip of the leading edge (regions U 
and III) by means of equation (A5) 

4. Estimation of the secondary corrections between the wing boundaries 
and, the reflected Mach lines from the trailing edge and tip, 
using as guides the trends of the primary corrections in steps 
2 and 3 

. Addition of the necessary lifting pressure in accordance with 
equation (l) 
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Where less accuracy is required, equation ( 7 ) and the quasi— conical terms 
in equations (11) and (A5) may he emitted. If the region between the wing 
boundaries and the reflected Mach lines is, small, then step 4 may also be 
excluded. 


Derivative for Damping in Pitch 


!Ehe derivative for the damping in steady pitch is calculated in the 
gginift manner as the pressure distribution; that is , the basic uncorrected 
value will first be determined based upon the expressions far the pressure 
distribution of tr iangular wings, and then primary and secondary correc- 
tive terms resulting from the cancellation of the excess pressure in the 
wake and outboard of the tips will be added. As before, the conical and 
quasi— conical components are developed separately in order to ascertain 
their relative magnitudes. 

If the distribution of pressure is known, as well aB the correspond- 
ing moment arm to the center of pressure of any sector of the wing, the 
damping moment in pitch far the basic plan form and for the separate 
regions affected by the cancellation of the exterior pressure can be 
readily calculated about some particular lateral axis. The spanwise axis 
passing through the center of gravity (assumed to lie an the x axis) haB 
been chosen In this analysis. However, since the pressure distribution is 
based on a pitching velocity about the y axis, the choice of any axis 
other - Hihti -the y axis requires a computation, of _ the additional angle— 
of-attack corrections given in equation (2). Hence, the damping deriva- 
tive calculated in the following sections is no;t, in itself, complete, 
except when the center of gravity lies at the leading— edge apex. 

Basic ywine far the damp'tng derivative .— With reference to figure 8, 
it is seen that an increment of damping moment due to the basic pressure, 
which is quasi conically distributed with respect to the y axis on an 
element of area (ds/da)da, is equal to the product of the resultant 
force dZ and the moment arm 5 e, the value of x being (—x^— x c# g # ) . 

Thus, from equation (3) and taking into account both halves of the wing 


dM = - 2 xdZ 


= _ 2 (la x — x V C g/ ga) gs f=aL. (W) da] (ll+) 

2 \4 A e.g.J |_ \J pv / 3 V 2V/ 0B c 6 ,/ m 2_a,2 Vda/ J 


In this equation, the resultant force is based on an average pressure 
equal to two— thirds of the pressure at the point x^,y^ at the end of the 
sector. The total moment is found by integrating this expression over the 
two ranges, 0 < a < a^, for which (from reference f) 
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_ = 

A mt-a 

dS _ m t 2c o 2 
da 2p(nn--a) 2 

and a.£ < a < m, for which 



dS _ ps 2 
da 2a 2 


Equatipn (l4) has been integrated numerically in Appendix B where it is 
expressed in the form of the damping— in-pitch derivative , that is. 



Since the pressure coefficient and the parameter ( — - ) are directly 

>> 2V-' 

related in the linearized potential theory, the derivative may be written 
as 



The tapered (m f ioj.) and the untapered (m.= m^) swept-back plan farms 
are treated separately. The result in each case represents the basic 
uncorrected value of the damping derivative, to which must be added the 
following primary and secondary corrections due to the cancellation of 
the basic and secondary pressure in the wake and outboard of the tips. 

« 

Primary corrections due to cancellation of basic pressure in wake.— 
If the trailing edges are subsonic, it is necessary to correct the basic 
value of the damping derivative for the effect of canceling the excess 
basic pressure in the wak9 by use of the oblique and symmetrical sectors 
of pressure. Considering first the portion of the basic pressure canceled 
by the symmetrical sectors,. the increment of force on an element of wing 
area (ds/dt 0 )dt 0 due to the cancellation of the conical component is 

dZi* =( |PV 2 ^ 5op(cp 0 ,ko) (J&-J dt D 
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That due to the quasi— conical component is 


*■' "l(K) F f 0 ^ a l{ r(<p ' 


k Q ) “ 


foJ"-o 


. 5d*>« Vo «„ 0 ,*o> i }&■) 
“tl V t 0 2 (l-^t 2 ) K 0 ^ 0 ° ° J V d W 


gin" 1 


dt. 


in which the average pressure is two— thirds of the maximum value at the tip 
(or leading— edge) end of the element of area. As illustrated in figure 9, 
the moment arms of the two respective forces are 

x i * = c o + g ^s -0 ©) — x c.g. 


xa* - c o + r (xs-^o) “ x c . g. 


giving the increments of induced moment 

dM x * = ?!» dZi* 


dMa’ = x 2 * dZ 2 * 

The total induced moment on the wing is found by integrating these two 
equations between the limits mj. and 1. If the rayB t Q intersect the 
leading edge as well as the tip, the integration must, be performed sepa- 
rately over the two ranges m^ < t 0 < hIq and ]% < w ^ ere “b# 

as shown in figure 1(b), designates TShe ray passing through the tip of 
the leading edge (i.e., mo=$ins/(Ps-mc 0 )). Far t Q < iUq 


x 8 



c 0 


dS = ps 2 
« Q 2t 0 2 


Com 

Xg = -2— + c r 


t Q -m 


dS _ 


n^Co 2 


dt Q 2P(t 0 -m)‘ 


while for t Q > m*-, 
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Based upon these values, the damping correction, expressed in derivative 
form, to he added to the basic derivative therefore becomes 


AC™ 1 = AC™ * + AC™ ‘ 


where 




(+s±\ r(i yisA) dt 0 + 

\ sy K L J V 3 Vo c Q y t 0 * Ko 

“t 

/ mc^Y n 1 / + 2 m _ fc ljSa . > ) 1 dt Q 1 (l 6) 

\Ps J J \ 3 ’t ; o“ nl C 0 / (■to - ®) 2 So J 


aln ~7§llS - } 


dt D + 


“o 


to 

®b 




dt r 


(17) 


If niQ is greater than 1, the second integrals in equations (16) and 
(17) are not required. The equations should be integrated by a graphical 
method. 

The anmp-ing correction due to the cancellation of the basic pressure 
in the wake in excess of that canceled by equations (16) and (17) is 
calculated in a s imil ar manne r. The induced m o ment on the wing due to a 
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single oblique canceling sector originating at a point A(x A ,y A ) on the 
trailing edge is first determined, and then the total moment is found by 
integrating the effects of all the oblique sectors superposed along the 
trailing edge. The two resultant forces, one (dZx**) arising from the 
conical portion of the canceling pressure and the other (dZ 2 **) from 
the quasi— conical portion, acting on an element of wing area (ds/dt)dt 
are shown in figure 10 far a single oblique sector. Considering first 
the case in which none of the rays from the oblique sectors intersect 
the leading edge, the average quasi— conical pressure is two— thirds the 
maximum value at the tip end of the element (for which x — x A is equal 
to P(s— y^/t), and the two increments of force from eqmtion (9) become 


“i" =GK) I (t^) 008 1 x " (i) 

az 2 " - ('ip/') da 

\2 / 3 It \ da / ,tx A mfc-a \ 


dt 


cos 1 X* 1 — 


t — 1— a 



dt 


The two respective moment arms are seen from figure 10 to be 


- _ .2 P(s-y A ) ' 

x i “ X A + 3 t ^-B- 


=, Xa + 1 g.^Aj _ x 

2 f A 4 t c, s* 

The term ds/dt is given in reference 7 as 


dS . 

dt 2a^ 2 t2 \m^.-a/ 


The two increments of induced moment due to a single oblique sector 
therefore become 


flMi" = xi»* dZx** 
dMa** = x 2 ” dZ 2 *» 


The total induced moment due to a single sector is now found by inte- 
grating these expreBBionB with respect to t over the range’ m^ < t < 1 
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on the wing. A second, integration with respect to a over the range 
-a^ < a < a^. then gives the total moment due to all the sectors super- 
posed along the trailing edges. Therefore the correction, expressed in 
derivative form, for the region adjacent to a subsonic trailing edge 
(regions I and III in fig. 1) becomes 


*V' m *Sx" + ^V’ 





x*» 



dtda + 



12 1 s ~y A m t~ t 

« 3 t x A n^-a 



t- 8 - »/ » » 2 

t-mt 1-a 



(18) 


•The results of the integration of equation (l8) with respect to t are 
given in Appendix C (case I); however, a graphical method of integration 
with respect to a is required for the complete solution. 

Configurations, in which the Mach lines from the canceling sectors 
in the lower range of a intersect the leading edge, require additional 
calculations for the range ny < t < 1, in which t = ny (ny < 1) is 
defined as the ray passing through the points A(x A ,y A ) and the tip of 
the leading edge (fig. 1(h)), that is. 


% 


= P 


8 -y A 

N _ x 
m X A 


Thus, the -integration with respect to t must he made over the two 
ranges m^ < t < ny and ny < t < 1. If only the conical pressure 
terms are considered, the area and center of pressure of an incremental 
sector in the range ny < t < 1 are 

dS _ xa 2 (hi — a ) a 
dt 2p(t - m) 2 


*L M 



+ i. m — a "N 

3 t - m J 


- x„ 
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This range disappears for those sectors which have Mach lines passing 
downstream from the tip of the leading edge, the limiting sector "being 
defined "by the ray (fig- 1(h) ) t where 

1 - Os/com) (1-m) 

a = 

L 1 - (Ps/m^mCo) (l-m) 

Hence far a > a. lt equation (18) applies. The damping correction far 
the region 0 < a < a^ is ! calculated also in Appendix C (case II), 

"but only for the c oni cal pressures since the quasi— conical pressures have 
"been shown to "be negligibly small. 

Thfl complete correction to the damping derivative resulting from 
the cancellation of all the "basic pressure in the wake is comprised of 
equations (16), (17), and (C2) (or (C6) or (C7), as required), in which 
the conical and quasi— conical terms are retained as separate components. 

PrlTrwry corrections due to can collation of "basic pressure outboard 

of tips .— Tha primary correction far a subsonic tip with zero rake can 
"be calculated in a mannar similar to that for a subsonic trailing edge. 
rpvm farces dZi*** and dZs 1 ** on an element of wing area at the tip 
due to the conical and quasi-conical components of the canceling pressure, 
together with the respective moment .arms, are shown in figure 11. The 
farce dZ 2 * 1 ’ corresponds to an average pressure equal to two-thirds of 
the Tnairi rmim pressure at the point of intersection of the ray t and the 
trailing edge, for which 


x - X/, 


s ( 1 _ jA 

t -mfc \ a a t/ 


From equation (12), 

az,". -(ip-v 2 ) 1 (^) 33003-1 X "’(S) " 

](i)« 


i 7 a v y i-x* » 

t(l+a) 


in which 


dS _ 

dt 2 
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\ a t V 


(™t - t ) J 
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The respective moment arms are 



Therefore, the primary correction to the basic damping derivative for 
regions H and HI at the tip becomes 


£C m »*« 

“q 


= * * * 
“qi 


£c m * * * 
“qa 




dtda + 



X*" - 


t-a 

t(l4a) 


Vi_x»«» 2 


. (S) 4t4a } 


(19) 


An analytical integration of equation (19) with respect to t in the 
range — 1 < t < 0 is given in Appendix D . A suitable graphical method 
for integrating the resultant expression with respect to a (as outlined 
in reference 7) is also included. The tapered and untapered cases are 
Considered separately. 

Secondary corrections As discussed previously with regard to the 
pressure distribution, the calculation of the damping derivative, to 
be complete, should include corrections for the secondary and smaller 
pressures. This is most simply accomplished by estimating the magnitude 
of the excess loading at the tip between the secondary (reflected) Mach 
lines a n d the edges of the wing (see diagrams of the pressure distribu- 
tion) and the distance from the pitching axis to the center of loading. 

For practical purposes, the excess pressure on the wing may be assumed 
to vary linearly from the maximum value along the wing boundary to zero 
along the Jfech line. The excess pressure along the tip is equal in magni- 
tude to the correction due to the cancellation of basic pressure in the 
wake; and that along the trailing edge, to the correction due to the 
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cancellation of "basic pressure outboard of the tips. The excess pressure 
upstream, from the leading edges may , in general, be neglected. This 
approximate method has been used in the following illustrative example. 

Illustrative example .- The derivative for the damping in pitch of 
the configuration shown in figure 1 has been calculated using equations 
( B 4), (C2), (D3), (16), and (17) and the procedure given in the pre- 
ceding p ar agraph for calculating the secondary correction. Far this 
example the center of gravity is assumed to be located at the apex of the 
l eading edge. The results are presented in table I, in which the conical 
and quasi— conical components are shown individually. 

It is apparent that the contribution to the total correction of the 
quasi— conical terms is only a small portion of the final value of the 
damping derivative and therefore may be dropped in many practical cases. 
The estimated secondary correction in some cases may also be insignifi- 
cant, particularly those in which the wing is considerably tapered. If 
such approximations are warranted, then possibly the simpler method of 
reference 6 for canceling the excess pressure in the wake nay be used. 

The calculation of the damping derivative therefore can be accomr- 
plished in the following steps: 

1 . Computation of the basic value for the over-all plan form by 

means of equation (B3) (or (b4)) 

2. Correction of the basic value for the effect of a subsonic 

trailing edge by means of equations (16), (17) , and (C2) 

(or (c6) or (C7), as required), dropping the quasi-canical 
terms where feasible 

3. Correction of the basic value for the effect of a subsonic tip 

by means of equation (D3), dropping the quasi-canical terms 
where feasible 

4. Estimation of the secondary correction for the excess secondary 

pressures as outlined, if required 

5. Addition of the angle— of -attack corrections in accordance with 

equation (2) 


This procedure follows closely those procedures given in reference 7 for 
the lifting case and in reference 9 for the rolling case. 


CONCLUDING REMARKS 


Through use of the method of the superposition of conical flows, the 
calculation of the pressure distribution and the derivative for the damping 
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in steady pitch, lias "been extended to include swept-back wings having all 
edges subsonic (provided they are also straight). The method presented 
is rigorous if carried out in detail, but for most practical purposes it 
can be shortened considerably by neglecting certain less significant 
terms in the final results. The method resulting from such simplifica- 
tion parallels closely those methods presented for swept-back wings in 
reference 7 for the case of lift and in reference 9 for the case of roll. 

Although this analysis does not include wings having positive or 
negative rake at the tips, it can be easily adapted to wings having nega- 
tive rake. The analysis is not complete for configurations in which the 
Mach lines originating at the trailing edge intersect the leading edge; 
however, the results given for such cases are believed to be sufficiently 
accurate for most applications. 

The analysis has shown that abrupt changes in pressure should occur 
along the various primary and reflected Mach lines an the wing, especially 
those originating at the tips of the leading edge. Such effects have 
been found previously for wings at constant angle of attack and in steady 
roll. It should be noted, however, that such discontinuities in pressure 
are not compatible with the flow that would be found on wing surfaces in 
a viscous, fluid. There is evidence that, because of the presence of a 
boundary layer on the wing surfaces, the pressure changes are much more 
gradual than those given by the foregoing theory. In particular, it is 
believed that a rigorous numerical analysis of the secondary corrections 
for the pressure and damping derivative (for region HI in fig. 1(a)), 
because of the moderating effect of the viscosity, would not be Justified 
from a practical standpoint. For this reason, approximate methods far 
calculating the secondary corrections were utilized in the present report. 


Ames Aeronautical laboratory, 

Watipnal Advisory Committee for Aeronautics, 
Moffett Field, Calif., May 12, 1950. 
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appendix a 


INTEGRATION OP PRESSURE CORRECTION DUE TO CANCELLATION OF 
BASIC PRESSURE OUTBOARD OF TIP 


The equation for the correction for the "basic pressure in the vici- 
nity of a streanwise tip clue to the cancellation of "basic pressure out- 
board of the tip (equation ( 13 )), namely. 





2m 2 -a 2 


( Py-ax) V (n^-a 2 ) (a-a Q ) ( 1+a) 


-da + 


5 gm2 ^ 2 da \ (Al) 

c o/ L a o( lm ) J (py-ax) */ (m 2 -a 2 )(a-a 0 )(l+a) J 
a o 

can he integrated in terms of elliptic integrals in the following Banner: 

Substituting — = gives 

c o c cP" 



--(SOi t i, W : *T 2i 


if 


m 2 


a(Py-ax)</(m 2 -a 2 ) (a-a Q ) (l+a) 
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m 
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fm+aHmr-a') 
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a(Py-ax) V (m 2 -a 2 )(ei-a 0 )(l4a) 



a ~ a o da + 

a o(l +a ) a(Py-ax)*/ (m?-a a )(l+ a )( a “ a o) 



a ~~ a o (min') (m-a) ^ 
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Equation 

tuting 


Thus, 
APp*' f = 


(A2) may be expressed in terms of elliptic integrals by substi- 


sn 
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where 
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The solution to equation (A^) reduces to the final form 



where 

* Q = [E(k) -K(k)] F(<P,k*) + K(k) E(q>,k«) 
% = tE(k) -K(k )3 E(+,k») + K(k) E(*,k») 
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The first "bracketed term in equation (A5) represents idle part of the 
correction due to the conical component of the sectors of canceli ng 
pressure; the second "bracketed term, that due to the quasi— conical com- 
ponent . 

Of particular interest are the corrections for points on the Mach 
line originating at the tip of the leading edge. Along this line the 
terms - a D and k in equation (A5) reduce, respectively, to m and zero; 
and, as a result, the second integral in equation (Al), corresponding to 
the quasi— conical component of pressure, vanishes (since a-a Q -3> 0 ). 

On the other hand, the conical component is finite along the Mach Utib 
and equation (A5) reduces to 



Thus equation (a 6) expresses an abrupt drop in pressure along the Mach 
line. 
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• APPENDIX B 

BASIC VALUE OF DERIVATIVE FOR DAMPING HT PITCH 


From equations (l4) and (15) the basic -value of the derivative for 
the damping in pitch may be written 
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After substitution of the respective expressions for x and dZ, equa- 
tion (Bl) reduces to 
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]} 


(B2) 


The final expressions for the two cases (in which > m) and 

after integration of equation (B2), are aB follows: 


\ 2 „2 


s -©O' c< 




(2m 2 +3j; 2 ) j/m 2 -a-fc 2 + 


3m 2 at 3 


(**<$ r JE3. pt*- 8 * 2 

\ ps / l Stm^^n 2 )^-^) L (m|.-a-|.) 2 2(m t 2-m2)(m fc -a t ) 


2nf t +12mPmjj 2 +m^ 4 1 m 3 ( 4m 4 — l Pm^ 2 m 2 + 23 m |. 4 ) ^ 

2(m t 2 -m 2 ) 2 -I 6m^. 3 (ulj- 2 -n? ) 3 


NACA TU 219? 


35 


m t m 2 (2m 2 +3ni fc 2 ) f _ x - _j_ 

— — COB — =— =■ — COS 

2(m fc a -«i 2 ) 3 */m t 2 -a a L m(mt-at) 



3 W°o /l \ p s / H 2 

/ni(.c 0 '\ 4 1” (2m 4 — 5mt 2 m 2 +in t Q a^+2m 2 m t a t ) A/m^-ot 2 

\ Ps / L 2(m t 2 -m 2 ) 2 (m t -a t ) 2 


m 3 (2m 2 — 5m^- 2 ) 3m-fc 2 ffl 2 

2m^ g (m^-rn 2 ) 2 2(mt 2 -Bi 2 ) 2 '^n4; 2 -ni 2 ^ m(mt-at) 


('cos -1 


cos 


■^)]} 


(B3) 


“ ° “t 


% = " (?) (fc) f {' 


f(2m 2 +a- fc 2 )v^-e^ 2 _ 
' 3^ 


\ Br / 


Ps ^ 105m? (m-a^.) 4 


j^( 64m 3 —46m 2 a^— 4®a^ 2 +a^_ 3 ) Vm 2 -at 2 — 64(mr-at) 4 J j* — 


4^ c . g .y /^o'n 

3 \ mc 0 / l \ Ps / 


*/m 2 -a^ 2 
Ps / a^. 2 


(?) - * («t) 3 ] } 


(b4) 


If the term x„ _ is zero in the above equations, then the angle-of— 

C*g* 

attack correction given by equation (2) is not required. 
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APPENDIX C 

CORRECTIONS FOR DAMPING DUE TO OBLIQUE . 
CANCELING SECTORS IN WAKE 


Thet corrections for the damping derivative due to cancellation of 
the basic pressure in the wake by means of the oblique sectors are given 
separately far the two cases >-l (i.e., Mach lines from trailing 

edge intersect tip) and < 1 (i.e., Mach lines from trailing edge 

intersect leading edge). 


Case I (n^ > l) 

" • / 

The correction for this case is given by equation (18), which, upon 
substituting the expressions for x* * and dS/dt, becomes 


' » * + * ‘ 

“q “qj. “qa 

Pm fc 2 s 2 t 


{£¥ / (¥9 OK/ h--- 

O o ®t 


2 P( S -?A^ 1 \ cos -1 X« »dt da + ^ P ^ * ■ f £_7A x 

3 t J t 2 3 Srtat 2 J ( mt _a) 3 x A \ da / 

th * ** 


& 




(Cl) 


Integrating equation* (Cl) with respect to t gives 

& ■• = (t) r H ai >A K-»> g f Aa-*c. s a _ 

""g da a(m-fc-a) 2 \ \ c Q / \ “t / 

• ! 0£) (^?) ♦ [ C 5 *^) + 
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(5 if { (*?-) [ - 

2m t -a(l+m t ) /(mfc-aHl-a) 1 ^ p / s-ja ^ |~ (Hi ( .-a)(2iii t +a) ^ 

ta^^l-a) v m^. J 8 V o 0 / L a 3 |!I 1 t 2 


mt-a 

2a 2 mj. 


amfc^Q+mfc 2 )-^ 2 / (m t .-a)(l-a) 1 \ 
lta 3 nij. 2 (l-a) v hLj. -* J 


(C2) 


where 


dPp 


A 


da 


/qc Q \ Qn^. T 2m 4 -3m a a 2 -ta 3 mj. 1 

' SB L (mt-a) 2 (nf-a 2 f /S ' J 


x = • 
A m-fc-a 


(C3) 


. JW_ 

A PCmt-a) 


Equation (C2) should he integrated graphically. Since the integrand of 
this equation is indeterminate for a = 0, the following expression ghnnid 
he used for the point a = 0: 


K "1 ■© ^ K(-Sf) 




Bsm 

2 1 ffic 0 m t 4 


{(^^•) ^mt(l-mt) + 


^ 6(l+mj.) (l-mt) 

c o _ 


3(l-^ t )(5+3mt) 


+ J^m^l+iiLj.) — 5(l+m fc ) 3 + '1 6 



«*) 
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Case II (n^ < 1) 


The correction for this case is composed of two parts, one for the 
range a 2 S. a 5: a t the °ther for the range 0 < a < a^. Thus, 



(*vO a> 



The first part for a > a^ is given by equation (C2). The second part 
has to he evaluated separately over the two ranges m. < t < m and 
“a < t < 1/ that is, o ^ 



If only the conical component of the canceling pressure iB considered, 
the correction far the range 0 < a < a j becomes 




{/T a [(— ) 

O 


+ 


I i cos ^ x ,,5a 


3 t 

' a 7 pi 


J 


2 2 

. 8 

2.2 




/ a t~ a \ 2 

VHrV 


ST 


“a 


( ^y co3 -^rr^(^f dtda} 
\ da / Jt 2p \ W . J 


dtda + 


(C5) 


When integrated with respect to t, this expression for a tapered wing 
reduces to 




\ 27 / 


fi. 

4 


\ a t / \mt-ay l\ C Q 



“Hr* , .1 ^%-a)(l-a ) coa ^L 2m t ^i a (l4 mt )' [ < g P(s-gr A ) [ % g ~a a _ coa -^L v , . *t g -* a , 

®t a «a v mt ma(l-mt) J 3 c Q ugan^* 2m t 2 a s 


2mt4a(l+mt) /(ifa)(l-a) ___-i 

i 1 - . i— ^VD 

4fla 2 m^ ^ 


210^-1% ( 1 -Hnt) , »J (l-a)(l-^n a )(m t -a)(m a Hm t ) 

T — ■ * ■ ■ . — _ 


1 

/ 



y iU M. J 


3*0 



r ffla— A 

— \ ( in n ^ rw_ mi n 

. ... . \ *»:—«/ i 

(m-a) 

L x A (m a -a)(m a ^a) n^-a J 


COB 


XI' _ 

1 % 


vhere 


yrr 


A. /I 

nir-a V (l-m)(mt-m) 


-/ 


a— o— 

*A 


^ «. « \ / 1 i — 

\m-<* 1 I XTlIH.-CiUJ 

+ " J'j _ n r— COB 


2(l-m)(mt-m 


r j°» v 


(ma-m) ( 1-m) (m^-m) J 


da 


X tt 
% 


V 


(l-a){ma-mt) - (»tfc-a) (l-i%) 
(1-^4;) (%-a) 


1 1 ~ \ /_ _ \ . \ / . \ 


(1-HBt) (%-m) 


(C 6) 


to 

S.O 
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For an untaper e& ving (11=10^), the correction la 


(s M ). 


fkLr ^ 

V zvJ 


§. (it ( !£±Y -f r^coa ” 1 X" - — 

4 \H/ \ ■ l \ c o / L “a® 11 xna 


1 /(m-a)(l-a) -i 2nwna(l-Hn) 

— - / — -■ gvu ■“ 

ita v m 


+ a p(^a) 
%(l-m) J 3 c 0 


v- a£ 

L 2 ion^ 2 a' 


V v f 

gup o - - 


m 2 -* 2 


a ‘ a * 2 m s a 2 


4 *a 2 m 


(m-a) 

3 *P 


'(£){ 


[ 6 + 2 ^. + 2 (^al 
L m fl -m l-m 


m 

- gwo 

^(l-m) 

2ttam 

J 

%-a 

r3(^A-^c. g .) 

2m*aMB a 1 

L (n^-m^m-a) 

1- X A ‘ 

®a“® -1 ■ 



- 

1 


tt _ 


^(!Ua)(l-i%) 'l 
(m-a)(iiVi-m) J 


da 
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APPEND IX D 

DAMPING CORRECTION DUE TO CANCELLATION 
OF BASIC PRESSURE OUTBOARD OF TIPS 


The correction far the damping ' derivative from equation (19) » after 

substituting the expressions for x^’”, ’ * * > and £±L and integrating 

with respect to t, becomes ^ 


AC_ » » * = AC„ * » * + AC„ ' * * 

m q “qx “qa 


. (t) «v « 

Q°o ^ 4 c n J a . \day\a a-fc / 


'o &i. 


( I _ EcjlELi?) -fjL 

V a p 8 / H 


]} 


A— + -A. 


/ a(l-Ki) ^ 

V mtCl+mt) L2mt 


1+2m t + 1 


( m t -a ) 2 “t-a v mt( 1+m t) '- 2 m t (l+m t ) m^-a 


M 

(5) 


da + 


P 2m t' 


s 




6 c f 


( 1 -fs-.s-M JL. 

\ a Ps / l 2m^ 


m-f-2a 

2(mt-a) J 


/ afl+a) 1 1 1+2m t /a(l+a) 

- A/mtCl+mt) -J Ml+mt)(m t -a) V “ 


4nL|.(l+mj.)(l+a) ** m^.(l+m^) 


/SiSLl + 3 _ l\ /_ _1 32^ 

v m^l+rn*)/ 4^'\ a ®t; / L 61^2 6 (mj--a ) 3 


/m^L r _ 

Vm+Cl+nu.) L 


l +2m t 


m-t+a+1 


_ + 3a-mt 

m^. ( 1+nij. ) L 8m^.(m^-a)(l+m^) 2 6(m t -a’) 2 (l+m t ) 6(m^-a) £ 


] 
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ks. 


1+2 ®t 


^4 (I+14) (l+a) V V 1+m t 


/ a(l+a) . „ I 

V 14(1+14) J 


da 


(Dl) 


M 


P 2 3 Ht £ 


(3) 


t [/ (^) 


2mt 

3 


fXit) j(a) to 


at 


] 


(D2) 


where 


J(a)= l(a-^) {(a~^) fe+ h(a) - 2 att (l^)(l4 a ) V^(^j'] 

j a(l+a) 1 1 

)v nu. ( 1+nu. ) J J 


ha_ jJn r _ i . . i+gm t 

4 \a a t / L 3mt 2 2m fc 2(l+m t ) 2 (l+a)'/ 114.(1+114-) 


and 

g(a) 

E(a) 

h(a) 

J(a) 


1 . f / a ,(l + a) _ ± 
mf-a |_V 14(1+14) 


14(1+14) 

1 L/ a \ . 1+2m t /a(l+a) ] 

14-a j_ s ^ 1 214(1+14) V 14 (1-H4) J. 

1+gm t j a(l+a) 1 

2(l-Hnt)*/mt(l+mt) J 

j a (l+a) 1 

+14)2 »/ i4(l-H4) J 


■ sjs |W* 

m^-e. [ J 


'a) g(a) - 




J(a) = ba-mt) g(a) + ?i2ii / a ^ 4a) 1 

3 (14-a) L 1+14 v 14(1+14). J 
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^3 


dPp 


UX Kn » 

Since the term — "becomes infinite as a approaches m, it is 

necessary to transform equation (D2), according to the procedure outlined 
in reference 7> to the following more suitable form for numerical integra- 
tion, where the two cases m p m^ (m^ > m) and m = m^ require separate 
solutions as noted: 


££„ ««* 
“q 


"(¥) 


!EV sL 

E- c o a 


■^I’fm) 2 m log^— — ^ _ 


I»(a) -S^m) .lnj 

a t 




/ 

*t 




/ 2 m 2 -a 2 > 

I da j- 

\ a > 

J 


(D3) 


Case I, m 3 ^ mj. 

E ’ (a) ■ (i - kJ {- ? fe + “ a >] + (i ■ - %“■) «■ w - *(i - p-w 


^t 

_2 


[^ + ' (a) ]}-?(r-4)[t + s(a) ](a- ? r t ) 


J* (a) = S ( 1 _ ic.gA L,( a ) 1 /a(l-fa) 1 1 _ 

2 \ a a t/ l\ a Ps J\_ Wtt t (l+m fc )(l4a)2y mt (l+m t )J J 

~t(~ ~ ~ ) [k(a) + 3 nU(a) + / a (l+a) 

2a a \a a+Z L ^ 2m, fl+m. } 2 f l 4 aV V m 


_3_ 

2 a 2 


2 m t (l+m t ) 2 (l+a)- v n^fl-mj.) J 

(i~$ { [£ +h(a) -- 


1 /-al-l.+a) ] 1 + ggt f 1 __ jV [" JJ(a) 

211^(1+11^) ( 1 +a) V m t (l+m t )Jj 8 \a a^ [ ^ 



kb 
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l +2m fc 


lJamt 2 (l+i ttfc ) 2 (l-m) 2 V 


/ a(l+a) 1 

V ^(l+m*) J 


in which g(a), p(a), h(a), ,)(a), and 2(a) have been previously defined, 
and 

gi(a) = —L— I" g(a) + — 2& " rl ■■ ■ 1 

2*/ am^(l+a)(l-HiL|.) -* 


,, * 1 r , , ,, . (l+2mt)(l+2a) 

p'(a) = p(a) + g* (a) + ; = 

“t -8 - L 4114. (I+114.) V am t (l+m t )(l+a) 


] 


h* (a) = — — 
m t~^- 


[(m-fc-Sa) 


g*(a) - 2g(a) + h(a) - 


(l+2mt)(l+2a) 


4 ( I+114. ) n/an4-( 1+a) (I+114.)" 


] 


J'(a) --i- j" J(a) + l>.(a) (l^m t)(naa) 

“t -3- L 8nL t (l+m fc ) ./am^l-fa) (I+114.) 


] 


71(a) = — i [32(a) + 3g(a) + (3aHm t )g*(a) + -i— /a(l+a) 

3(mt-a) L I+114. ,7114.(1+114.) 


(l+2a)(mtta+l) 


2(1+114;) J am t (l+m fc )(l+a) 


] 


Case H, m = 04. 


E ’ (m) - (l-^) {-s?[i + s(m) ] + (s-Ti^) B ' <m) -.f(i-^) p ’ (m) 

4 

;[i + l(m) ] }-^(i-H)[i +g(m) ] (i"V) 
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J* (m) 


where 
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(?■£) [ hw + + ] - 


K-"^) [ J ’ (m) + dfe ] 


g(m) = 


2m+l 

2 m(l+m) 


g*(m) 


1 

8m a (l+m) 2 


p(m) = - 3+8sB-8 m 2 
8 m 2 (l+m ) 2 


p*(m) = -| +g * „ 

8 m 3 (l+m ) 3 


h(m) = — — + 3. 

m 8 m(l+m) a 


h'(m) = i 

Sm^fl+m ) 8 


j(m) = 4+15m+2lta a +8m 3 
2 l|m 2 (l+in ) 8 


3 '(m) 


3 + 8 m 

64 m 3 (i+m ) 4 


Km) = -5*i_ 
^(l+m ) 2 


j»( m ) _ _ l+3gH-12m a +8ni 3 
24 m 2 (l+m ) 3 
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TABLE I.- CALCULATED DAMPING DERIVATIVE FOR WING SHOWN IN FIGURE 1 


Quantity 

Magnitude 

Percent 
of total 

C_ (uncorrected) 

-12.174 . 

179.3 

AC m 1 * (T.E., conical) 
^li 

1.063 

-15.7 

AC_ 11 (T.E., quasi conical) 
^12 

.046 

-7. 

AC™ * (symmetrical, conical) 

Ml 

1.883 

-27.6 

AC™ * (symmetrical, quasi conical) 

Iu q2 

.129 

“1*9 

AC™ 1 ** (tip, conical) 
^li 

2.431 * 

-35.8 

AC ,,T (tip, quasi conical) 

\l2 

.231 

-3.4 

Estimated secondary corrections 

— 4 oo 

5.9 

Cmq (corrected) 

-6.791 

' 100.0 









MCA TN 2197 


k9 



Figure /. - Wing plan form , coordinate system, and prin 
cipai symbols used in analysis. 





(b) Coordinate system and principal symbols 


Figure /. - Concluded. 
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(a) Pressure field canceled by symmetrical sector. 





Figure 2. - Field of pressure in the wake to be canceled. 
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Figure 3.- Symmetrical pressure canceling sector in the wake . 
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O 20 40 60 80 100 


Percent chord 
(a) Section A- A. 

Figure 7- Pressure distributions along various sections of the illustra- 
tive plan form in steady pitch. 
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_2 Quasi-conica! correction (TE)' 
Quasi-conicai correction (tip)- 
-. 4 - .Conical correction (TE)- 
Symmetricai conical correction- 
-q. .Conical correction { tip )- 


Estimated secondary corrections ( TE and tip) 


-8 


20 


40 60 

Percent chord 

(b) Sec tion B-B. 


Figure 7.- Continued. 
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Figure 9. — Forces on an element of wing area due to the conical and quasi- 
conical components of the symmetrical pressure canceling sector in Hie 
wake. 
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Figure !0- Forces on an element of vting area due to the conical and 
quasi- conical components of an oblique pressure canceling sector 
in the wake. 
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Figure //. - Forces on an element of wing area due to the conical and 

Quasi -conical components of a pressure cance/ina sector at the tio. 
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